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We propose the assumption of quantum mechanics on a discrete space and 
. time, which implies the modification of mathematical expressions for some 

postulates of quantum mechanics. In particular we have a Hilbert space 
■ where the vectors are complex functions of discrete variable. 

^ . As a concrete example we develop a discrete analog of the one-dimensio- 

nal quantum harmonic oscillator, using the dependence of the Wigner func- 
tions in terms of Kravchuk polynomials. In this model the position operator 
has a discrete spectrum given by one index of the Wigner functions, in the 
£h ' same way that the energy eigenvalues are given by the other matricial index. 

Similar picture can be made for other models where the differential 
q-i equation and their solutions correspond to the continuous limit of some 

difference operator and orthogonal polynomial of discrete variable. 



Oh 



X 

C<3 1 1- INTRODUCTION 



In the last Symposium on Fundamental Problems in Quantum Physics 
[1] I explored the hypothesis of a realistic interpretation of lattice theories 
based on some ontological model that presupposes some fundamental net- 
work previous to the concept of space and time. According to this model 
the structure of space-time is a consequence of the relations among these 
fundamental entities, and gives raise to a discrete character of the space and 
time variables. 

Even the standard interpretation of quantum mechanics in our model 
is conserved the assumption of a discrete space and time introduces some 
drastic changes in the mathematical formulation of quantum mechanics. 
These consequences have been also shared by recent authors who use lattice 
models as a mathematical tool. In particular non-commutative geometry 
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leads in a natural way to difference operators [2]. Quantum groups and q- 
analisis implies a deformation of space-time groups with a q-casimir denned 
on a space or time lattice [3]. Lattice field theories are widely spread in an 
attempt to overcome infinities in perturbation theories [4] . Special functions 
and orthogonal polynomials on continuous variable are studied together with 
those of discrete variable [5] . Non-euclidean crystallography requires the use 
of discrete groups of hyperbolie type that has been recently developped [6]. 
Recent literature advocates for the application of discrete models to unify 
general relativity and quantum mechanics as Wheeler, Ponzano and Regge 
and others have proposed [7]. For this purpose some modern tools are used, 
such as discrete topology and partition functions defined over symplicial 
networks. 

2. THE POSTULATES OF QUANTUM MECHANICS ON Z 4 

The assumption of discrete space and time imposses some restrictions 
on the mathematical expressions for the postulates of Quantum Mechan- 
ics. The question is now whether this assumption keeps the analogy with 
standard formulation in the continuum case (in the limit both formulations 
coincide) and at the same time avoids the unwanted infinities. The answer 
is, generally speaking, in the affirmative. 

In order to be more explicit, we start with the Hilbert space: this must 
be defined over the gaussian numbers (complex of integer components) with 
scalar products constructed with summation instead of integrals 



where we use the notation fj = f (je), e being the fundamental length of 
the one-dimensional lattice, j integer numbers. 

As an example, we take the following orthonormal basis 



b 



j=a 





N - 1 



where 



2 Trm 



m = 0, 1, • • • N - 1 



m 




satisfying 



N-l 




j = o 
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with respect to which a finite Fourier transform (equivalent to a Fourier 
series) can be defined 

N-l N-l 

Fj = ^ ' a mfj (km) ! O'va = ^ ' fj (^m) Fj 
m = J = 

Notice that although the space-time variables are discrete the functions 
can be continuous. 

If we consider observables, the corresponding operators must be self ad- 
joint with respect to the scalar product mentioned before, and the spectrum 
is always discrete. As an example, we give the eigenvalues and eigenfunc- 
tions of the position and momentum operators. 

X : jS jt = £S je , I fix ; j = 0, 1, ■ ■ ■ N - 1 



P : i-Ajfj (k m ) — - — fj (k m ) 



where 5ji is the Kronecker function and fj (k m ) is defined as before. 

With the help of the scalar product we can defined espectation values, 
projection operators, density matrix, uncertainties or mean values of some 
operators. 

Suppose now that a physical system is represented by a state vector ip (t) 
depending on discrete time (t = nr). If H is the operator corresponding to 
the hamiltonian of the system (for simplicity we take H constant in time) 
the Schrodinger equation for that system is 



1 _ 

T 



the solution of which is 



^ = (irfSf ) ^ 

with the initial condition ipQ. 

Here A n is the forward difference operator A n ip n = ip n +i — tp n and A n 
the mean operator A n tp n = | (tpn+i + VVi) 

As in the continuous case we can define an unitary evolution operator 

which is unitary (because H is self adjoint) and satisfies the difference equa- 
tion 

-A n U n = HA n U n 
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If we use the Heisenberg picture the evolution of some operators in time 
is given by 

A n = U+A U n (1) 

where A n is some operators depending on discrete time and U n is the evo- 
lution operator defined before. The Heisenberg equation now reads [8] 



- AA n — - 1 : [A n , H n ] i . 

1 — 2 ZT H 1 + 2 lT ^ 



the solution of which is (1). 

An other scheme may be used if we take the symmetric difference op- 
erator 



S n A n = (A n+1 _-A n _i), 



then 



i 1 



Also we have 

i 1 ( t 2 \ 1 

r {An+i ~ An ~ i) = i+ir^ 2 K; Hn] 1 + 

Some simplification can be achieved if we take the particular case H 2 = 
[t can be easily proved the following equation for the operator in the 



1. It can 

Heisenberg picture: 



1 AnA n = ^^( 1 -^)[A n ,H] 



Lorente 



5 



In the last three equations the dependence on the Hamiltonian operator 
H is linear as in the continuous case. 

The realization of operators in the coordinate or position representation 
is given by the substitution 

The substitution is not unique. We will discuss in the next section 
some different realization for the position and momentum operator of he 
harmonic oscillator. 



3. QUANTUM HARMONIC OSCILLATOR OF DISCRETE 
VARIABLE 

The quantum harmonic oscillator is described by the Schrodinger equa- 
tion 



~2~ 



with^^A t = as , a=y/*fr, A = jg 
For simplicity, we take a = 1 
The normalized solutions are 

Y> n (s) = (7rh n n\)~K- a2 / 2 H n (s) , n = 0,1,2,- ■■ (2) 

where H n (s) are the Hermite polynomials. 

The ip n (s) are eigenfunctions corresponding to the eigenvalues A = 
2n + 1, and they satisfy the following recurrence relations: 



i) 2s^J n (s) = ^2 (n+1) Vn+l (s) + V 7 ^ («) (3) 

ii) («) = — \/2 (n + 1) V'n+i (s) + V2n V» n -i (a) (4) 

From these two relations one defines the creation and annihilation op- 
erators: 



a f Y>n (s) = -7= f S - ) 0) = + l-0n+l («) (5) 



\/2 V ds, 

dlpn (s) = (s + («) = vWn-1 (s) (6) 
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It is well known that the Hermite polynomials are the continuous limit 
of the Kravchuk polynomials of discrete variables k n (x) and the weight 
function of the Hermite polynomials is the continuous limit of the binomial 
distribution p (x) which in turns is the weight function of the Kravchuk poly- 
nomial [9] . But the product of the Kravchuk polynomials time their weight 
function is proportional to the Wigner functions d 3 mm , (0) that appear in 
the generalized spherical functions 

D lnm> («' Z 3 ' 7) = exp (-ima) d 3 mm , (0) exp (-im'j) , (7) 

namely, 



dim' (P) = (-ir- m '^p^) (x,N) (8) 

where d n is some normalization constant, and n = j — m , x = j —m' , p = 
sin 2 ((3/2). 

This connection betwen the functions of discrete and continuous vari- 
able suggests that the solution of the quantum harmonic oscillator are the 
continuous limit, up to a factor, of the Wigner functions. In order to prove 
this Ansatz we compare the recurrence relations of the two types of func- 
tions as it is done for the orthogonal polynomials of discrete variable. In 
our case we take the differential equation for the Wigner function [10] 

±± d j (0) + m '~ mcos/3 <$ (3) - 



= y/UTm) (j±m + l) d j m±1 , m , (J3) (9)(10) 

From this we deduce two recurrence relationes: 

., n m' — mcos ,j , „. 
1) ^ sin/? P dl m , (0) = 



= yjti-m) (j + rn + l)d? m+lim , (0) + 



+VC? +m)(j-m + l)dl_ hm , (0) (11) 



ii) VC? + m ') 0' " m ' + l ) d m,m'-l (P) ~ 



-y/ij-m!) (j + m' + l)< m , +1 (P) = 



= yj (j -m){j + m+ l)d ] m+l m , - 

-^{j + m) (j- m + l)dl_ l m , (12) 
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Note that the last expression has been obtained with the help of the 
well known property of Wigner functions 

We suppose that 

^lim^ C n (N)di m ,(P) = i> n (s) (13) 

where we take m = j — n, m'=j — x, N = 2j, x = Np + \J2Npq s, and 
C n (iV) some normalization constant to be determined. 

We compare the recurrence relations i) that is to say, formulas (3) and 
(11). We divide the second one by yfj and substitute d 3 mm , {(3) by , 

withv n (x) = C n (N)di_ nd _ x (f3). 
The result is: 



j — x — (j — n) cos (3 v n (x) / 2n (N — n + 1) t>„-i (x) 

VJsin p C n {N) ~ V 2j Cn-i (x) ' 

\ 2{N-n) (n + 1) v n+1 (a) 
V 2j C n+1 (x) 



( (2p — l)n\ w„ (x) / : ~ / rT v„+i (x) 



N C n (N) 

In the limit N — > oo this expression goes to the recurrence relation (3) 

P™ vided S§y = cS§y = X > or Cn W = const = L 

The recurrence relations ii) formules (4) and (12) can be compared by 
the same method. We substitute 

v n {x) = di_ n>j _ x {(3) (14) 



in (12) and divide both sides by y/N. 
The result is 



( x + i) Vn (x + 1) - Wx^ V (x - 1) 
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Substituting x = Np + ^2Npqs, and extracting yJ2Np = ^ in the left 
side, we obtain 



N 



1 + \/|; s+ ^)"» (x+1) - 



= J2n (l - !L j^)»»-l to " W 2(n + 1) (l - (i) 

In the limit N — > oo, /i — ► this expression goes to 

lim - {V'n (s + /i) - (s - /i)} = -J2(n + l)tp n (s) + v/2n^n-i (s) 
h —* fl v 

that coincides with (4) 

We can use these results to construct creation and annihilation opera- 
tors for the Wigner functions. We define 



1 \ m'-mcos(3 



4j 



(j - m') (j + m' + 1,^ 
\j-m) (j + m + 1) j 



4i 



2j 



m+l,m' 



1 ( m' — mcos (3 j 
V2{ v(?sin dmm ' 



(j + m') (j — m' + 1) 
4j 



+ w (y '- TO/)( ?' + m/ + 1) < H , +1 (^) = 



4j 



(j + m) (j — m + 1) 



2j 



<4-l.m' 09) 



(15) 



(16) 



Using the limit of the recurrence relations we obtain 
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N — > oo 





i s + Ts 













(17) 



^lim^ AUl m , (/J) = _ ( s - _ ) ^ ( a ) = aHn (s) (18) 

Relations (15) and (16) suggest that the creation and annihilation oper- 
ators are connected with the raising and lowering operators for the spherical 
harmonics Yj m . In fact, multiplying (15) and (16) by Yj m i and adding for 
m' we have 



m 

or 



AY jm = ^j\Jti- m) (j + m + l)Y j>m+1 = -±=J + Yj m (19) 

similarly, 



or 



A*Y jm = -L=yJ(j + m ) U~m + l)Y jt m _i = ~^=J-Y jm (20) 

In order to make more transparent the connection between the creation 
and annihilation operators with the raising and lowering operators of the 
spherical harmonics, we take the commutation and anticommutation rela- 
tions of the former operators. 



[AA^-^A)Y jm = ^(J+J--J-J + )Y j 



2j 

1 m / n \ 

-—2J z Yj m = —Yj m = I 1 — — 1 Yi'm 
2j J V J , 



Substituting Y, m = £ c^ m , (/3) Y^/ we get 

m' 

AAt]^ mm ,(/3)=(l-^^ m ,(/3) (21) 



which in the limit j — ► oo goes to 
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Ipn (s) = (s) 



Similarly 
(aA* + A+a) Yj. 



or 

(aA* + A*a) d? mm , (/?) = | (2n + 1) - ^ | < m , (/?) (22) 
which in the limit j — ► oo goes to 

(act* + a* a) Vn («) = (2n + 1) tp n (s) 

If we multiply both sides by hu/2 we obtain the eigenvalue equation for 
the hamiltonian. 

The interpretation of this model can be taken from the quantum har- 
monic oscillator of continuous variable. The energy levels are equally distant 
by the amount Tioo and are labelled by n = 0, 1, 2, • • • oo. In the quantum har- 
monic oscillator of discrete variable we have also the discrete eigenvalues of 
the hamiltonian connected with the index m = j — n of the Wigner function 
d 3 mm , (P). These values are equally separeted but finite (m = — j, • • • + j). 

Similarly the eigenvalue of the position operator A+A + are also discrete 
and connected to the index m' = j — x of the Wigner functions but finite 
(W = 

The integer numbers x = 0, 1, • • • 2j are related to the quantity x = 
as where s is the continuous variable and a = y/Mu/fi. Since x is a 
pure number and s has the dimension of a length, the spacing of the one- 
dimensional lattice is equal to 1/a = \JhjMui. Therefore the Planck's 
constant H play an role with respect to discrete space similar to the role 
with respect to discrete energy values. 



a, a 



— (J+J- + J-J+) Yj rn — - (j 2 — J%\ Yj m — 

^3 3 

j (j (3 + 1) " rn 2 ) Y jm = j(2n + 1) - j j Y 3 



24. CONCLUDING REMARKS 

The analysis we have made for the quantum harmonic oscillator of dis- 
crete variable can be applied to other cases, where the functions involved 
are orthogonal polynomials of continuous variable the limit of which are 
some orthogonal polynomial of discrete variables; we give some examples: 
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1. The function fj (k m ) described in section 2 are polynomials of dis- 
crete variable the continuous limit of which are the exponential function. We 
have developed a new scheme for the Klein-Gordon and Dirac field equation 
that can be extended to lattice gauge theories [11]. 

2. The solution of the Schrodinger equation for the hidrogen atom is 
given in term of the orthogonal Laguerre polynomials and the spherical har- 
monics. The radical equation can be translated into the difference equation 
for the Meixner polynomial of discrete variable. 

3. The quantification of the electromagnetic fiels leads to the D'Alam- 
bert equation the solution of which are given in terms of the Bessel spherical 
functions that are related to the trigonometric functions. These functions 
suggest the connection with the orthogonal polynomials of discrete variable, 
that are solutions of difference equations of the hypergeometric type. Gen- 
eral speaking a parallel study of discrete and continuous model can be made 
similar to that made by the russian school of mathematicians with respect 
to orthogonal polynomials. 
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